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A simulation approach is developed to obtain the linear energy transfer (LET) spectrum of all secondary 
ions and predict single event upset (SEU) occurrence induced by neutron in memory devices. Neutron reaction 
channels, secondary ion species and energy ranges, and LET calculation method are introduced respectively. Ex- 
perimental results of neutron induced SEU effects on static random access memory (SRAM) and programmable 
read only memory (EEPROM) are presented to confirm the validity of the simulation results. 
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I. INTRODUCTION 


Single event upset (SEU) induced by neutrons in advanced 
memory devices is currently recognized as one of the major 
concerns in issues of reliability [1—4]. Although neutrons are 
uncharged particles and can’t directly induce electron-hole 
pairs in semiconductor devices, they transfer energy and pro- 
duce recoil and spallation products that are highly ionizing 
and therefore capable of creating electron-hole pairs in semi- 
conductor materials [5, 6]. Secondary ions, such as recoil 
and spallation products, can create large numbers of excess 
carriers and possibly cause SEU in memory devices. Pre- 
vious literature reports the phenomenon and mechanism of 
SEU [7, 8]. Experiments show that neutron induced SEU has 
been identified in memory devices and the cross sections for 
various memory devices are measured in different neutron 
radiation environments [9]. Simulation of neutron induced 
SEU effects is mainly based on the Monte Carlo method, as 
this method is useful for simulating SEU of a memory chip 
induced by neutrons and providing a detailed description of 
the random processes of SEU. Li H’s paper provides a Monte 
Carlo method to simulate the process of energy deposition in 
sensitive volumes and to calculate the SEU cross-section [10]. 
In addition, as one of the Monte Carlo toolkits, Geant4, is 
currently playing a more and more important role in many 
nuclear technology applications, and is especially used for 
neutron transport and nuclear reaction simulation. In Pete Tr- 
uscott’s paper, Geant4 is performed to simulate neutron SEU 
effects. SEU rate predictions are in good agreement with 
experimental data for large (~ 300 um) to small (~ 0.5 um) 
feature-size devices [11]. All the above-mentioned research 
provides a universal understanding of neutron induced SEU 
effects, and demonstrates that the Monte Carlo method, such 
as Geant4, is a reasonable and effective approach for neutron 
SEU simulation. 
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However, few research papers focus on how to predict the 
occurrence of SEU for different kinds of memory devices in 
neutron radiation fields. Practical neutron radiation, such as 
the reactor neutron, terrestrial neutron, and atmosphere neu- 
tron, provides both thermal and fast neutrons over a wide 
range of energies. In this paper, the experiments are per- 
formed at the Xi’an pulsed reactor (XAPR) at the state key 
laboratory of intense pulsed radiation simulation and effect 
in China, so in the simulation the incident neutron spectrum 
is the XAPR neutron spectrum. Because of the complicated 
device structure and neutron spectrum, it is quite challenging 
to determine whether memory devices exhibit SEU effect in 
the neutron radiation environment. 

It is a well known fact that linear energy transfer (LET), 
referring to energy loss of an incident ion per unit length in 
material, is an important parameter in detailed analysis of the 
SEU effect. It is commonly presumed that SEU occurs when- 
ever the LET of ions is more than the SEU threshold of the 
memory device. According to this knowledge, it is of the 
essence to calculate the LET range of all secondary ions in- 
duced by the neutron. In this paper, the objective is to predict 
whether SEU occurs in certain memory devices exposed to 
neutron radiation. First of all, in order to determine the en- 
ergy spectra of each secondary ion induced by the neutron, the 
Geant4 Monte-Carlo simulation program is adopted to pro- 
cess the transport of neutrons in silicon dioxide and simulate 
the interaction between incident neutrons and target atoms. 
Secondly, the LET of secondary ions is calculated using Stop- 
ping and Ranges of Ions in Matter (SRIM). Thus SEU will 
occur when the maximum LET is more than the SEU thresh- 
old of a memory device, and vice versa. Finally, reactor neu- 
tron radiation experiments on static random access memory 
(SRAM) and programmable read only memory (EEPROM) 
are performed to confirm the simulation prediction. 


Il. SIMULATION PROCEDURE 


In our simulation, in order to predict the occurrence of 
SEU induced by neutrons in silicon-based memory devices, 
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the problem is to calculate the LET spectrum of all secondary 
ions, considering the incident neutron spectrum, target ma- 
terial, and nuclear reaction channels. The incident neutron 
spectrum defines the primary neutron energy in the simula- 
tion by Monte Carlo random sample. The target material is 
defined as silicon dioxide. The top layer materials, such as 
copper, are not considered in the simulation for two reasons. 
The first is that in silicon-based devices, silicon dioxide is the 
major material for insulation and isolation of adjacent tran- 
sistors, while other layer materials are minor. The second 
is that the sensitive volume material is silicon, so the reac- 
tions between the sensitive volume material and neutron can 
be included in the reactions between silicon dioxide and the 
neutron. For the reactions between the neutron and silicon 
dioxide, there are 30 comprehensive nuclear reaction chan- 
nels, as shown in Table 1. 


For analysis presented within this paper, the prediction 
simulation is primarily based on the Geant4 and SRIM sim- 
ulation toolkits. Type and energy distribution of secondary 
ions, which are exported from the Geant4 simulation, are im- 
ported into SRIM to calculate the LET value of each sec- 
ondary ion. 


Geant4 uses the NeutronHP model for neutron interactions 
from thermal energies up to 20 MeV. In this paper, the Geant4 
Monte-Carlo program is first adopted to deal with collision, 
moderation, and nuclear reaction, and simulate the energy 
distribution of secondary ions produced from the nuclear re- 
action. The initial neutron energy spectrum, physical pro- 
cess, and reaction channels should be taken into considera- 
tion. Figure | shows the initial neutron spectrum provided 
by XAPR. In the Geant4 primary generator action file, the 
primary energy of the incident neutron is sampled from the 
spectrum mentioned above by random sampling. The neu- 
tron energy ranges from 3.6 x 107? eV to 20 MeV. The ma- 
jority of neutrons are in the low energy range, and differential 
flux decreases with an increase in neutron energy after 10 eV. 
In Geant4, neutron transport at energies below 20 MeV is 
based on empirical elastic and inelastic neutron cross-section 
data drawn from many of the standard database sources [12]. 
Ionization due to the recoil nucleus and secondary charged 
nucleons includes accurate treatment of low-energy electro- 
magnetic processes, using electronic stopping power for pro- 
tons and a particles, multiple scattering, and 5-ray production 
from ions [13]. 

In silicon based semiconductor devices, when the initial 
neutron is passing through a sensitive region, the neutron in- 
teracts with nuclei from the device material, predominantly 
28Si and !6O. The reaction is elastic when conserving mo- 
mentum, or inelastic when the neutron is absorbed followed 
by an ejection of nuclear fragments. The characteristics (ion 
type and energy) of the secondary ions responsible for SEU 
strongly depend on the type of nuclear reaction [14]. As 
shown in Table 1, the various nuclear reactions can only oc- 
cur above a certain energy threshold of a few MeV. As neu- 
tron energy increases, the number of reaction channels in- 
creases. Consequently, neutrons with low energy will not be 
very efficient in creating SEU effects through limited reac- 
tions. For neutrons above 17.8 MeV, all of the reaction chan- 
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TABLE 1. Schematic of neutron reaction channels and associated 
threshold energies 


Serial number Reaction channels Neutron threshold energy 


1 >8Si(n,2n)’ Si n> 17.8MeV 
2 °8Si(n,na)“*Mg n> 10.4MeV 
3 8Si(n,n) "Si n> 1.8MeV 
4 >8Si(n,g)’Si 0 

5 >8Si(n,p)°Al n> 4MeV 

6 *8Si(n,d)?’Al n> 10MeV 

7 >8Si(n,a) Mg n> 2.7 MeV 
8 °Si(n,2n)*Si n> 8.8MeV 
9 2Si(n,na) Mg n> 11.5 MeV 
10 2Si(n,np) Al n> 12.3 MeV 
11 °Si(n,n) Si n> 1.3 MeV 
12 *°Si(n,g)*?Si 0 

13 *°Si(n,p)”’Al n> 3MeV 

14 Sin, a) Mg n> 1.3 MeV 
15 30Si(n,2n)”Si n> 11 MeV 
16 30Si(n,na) Mg n> 12MeV 
17 30Si(n,np)” A1 n> 14MeV 
18 30Si(n,n) Si n> 2.3 MeV 
19 30Si(n,g) Si n> 1.5 MeV 
20 30Si(n,p)® A1 n> 8MeV 

21 30Si(n,a) Mg n> 4MeV 
22 °O(n,2n)5O n> 17MeV 
23 °O(n,na)?C n> 8.5 MeV 
24 °O(n,np) ÎN n> 13MeV 
25 ĉO(n,n)!™O n> 6.5 MeV 
26 ĉO(n,g) "O 0 

27 °O(n,p) SN n> 10 MeV 
28 °O(n,a) 2C n> 4MeV 

29 °O(n,d) SN n> 10 MeV 
30 O(n,t)'4N n> 17MeV 


$ 


Differential Neutron Flux (n/eV/cm*/s) 


Neutron Energy (eV) 


Fig. 1. Reactor neutron spectrum at XAPR. 


nels are open, increasing the probability of generating a sec- 
ondary ion able to induce a bit upset. 


The Monte Carlo simulation code named SRIM (formerly 
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TRIM) is a collection of software packages developed by 
Ziegler and Biersack, which is written to permit the calcu- 
lation of ion deposition profiles in materials exposed to ener- 
getic beams of ions [15, 16]. SRIM tables of LET, stopping 
power, and projected range versus ion energy are especially 
useful for computation of ion transport in a wide range of sim- 
ulation applications. In this paper, SRIM is used to calculate 
the LET value of secondary ions based on ion stopping and 
range tables. In the LET calculation, atomic number, mass 
number and energy of the incident ion, and type of target ma- 
terial should be defined. Figure 2 shows the LET of eight 
types of ions, which are secondary ions of neutron reactions 
with Si and O atoms. In silicon, the energy to create a pair of 
electron-holes is 3.6 eV. Therefore, 1 MeV cm?/mg LET pro- 
duces about 10.3 fC/um of silicon. As shown in the figure, 
the LET of each ion increases as the energy increases initially 
and decreases afterwards, and the ion with a larger atomic 
number can evidently induce a larger LET. 


150 
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a 
Oo 


lon Energy (MeV) 


Fig. 2. (Color online) LET versus energy for eight types of sec- 
ondary ions. These plots were generated using SRIM. 


The simulation results of the energy range and maximum 
LET are shown in Table 2. The energy range of each sec- 
ondary ion is simulated by Geant4, and maximum LET is 
calculated by SRIM. The maximum LET of neutron induced 
secondary ions is 9.28 MeVcm?/mg, which is caused by 
Mg at 3.79MeV. It should be noted that the maximum 
LET for protons and « particles in silicon is, respectively, 
0.54 and 1.45 MeV cm?/mg. In addition, the SEU cross sec- 
tion of these particles is smaller than most of the other sec- 
ondary ions. By dividing the LET range into successive in- 
tervals and calculating the proportion of ion numbers in ev- 
ery LET interval, we can obtain the LET spectrum of neu- 
tron induced secondary ions, as shown in Fig. 3. The LET is 
mostly distributed from 0 to 3 MeV cm?/mg. The peak point 
at 0.4MeV cm?/mg is mainly caused by secondary low en- 
ergy '°O, generated from the reaction of the primary neutron 
with material oxygen atoms. The LET spectrum is highly 
useful for predicting the SEU occurrence and further investi- 
gating ionizing effect on the electrical response of the mem- 
ory cell. In this paper, LET simulation results are used to 
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TABLE 2. The energy and maximum LET of secondary ions in- 
duced by reactor neutron 


Ion species Energy (keV) Maximum LET (MeV cm*/mg) 


Si 0-2080 8.16 
Si 0-810 4.49 
305i 0-1360 6.19 
3193 0-10 2.20 
7A] 0-2130 7.01 
8 AI 0-1930 6.52 
30A] 0-870 3.83 
Mg 0-3790 9.28 
25Mg 0-3840 9.18 
Mg 0-3460 8.70 
Mg 0-2460 7.31 
160 0-3190 7.04 
9 0-1510 5.61 
My 0-10 1.03 
ISN 0-2820 5.98 
16N 0-2310 5.81 
2G 0-4990 5.14 
BC 0-5000 5.14 


judge whether SEU will occur by comparing the maximum 
LET and the threshold of the memory device. 


Normalized Distribution of Secondary lons 


LET (MeV-cm°/mg) 


Fig. 3. LET spectrum of neutron induced secondary ions, which can 
be used to estimate SEU rates. 


Previous literature shows that the SEU threshold is 0.3- 
0.7 MeV cm?/mg for CMOS SRAM with a feature size of 
0.18-0.5 um, while the maximum LET is calculated to be 
9.28 MeV cm?/mg in the reactor neutron radiation environ- 
ment [17, 18]. For this reason, it can be confirmed that SEU 
will occur in SRAM exposed to reactor neutron radiation. 
Previous research shows that the SEU threshold is approxi- 
mately 37 MeV cm2/mg for EEPROM 28LV010, and is much 
larger than the maximum LET induced by the neutron [19]. 
Therefore, it can be confirmed that the reactor neutron will 
induce no SEU in the EEPROM devices with similar techno- 
logical processes. 
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I. EXPERIMENTAL PROCEDURE 


Neutron SEU experiments were performed at XAPR. The 
nuclear reactor facility is located at the Northwest Institute 
of Nuclear Technology (NINT) in China, and provides both 
thermal and fast neutron beams over a range of energies. Neu- 
tron fluence was approximately 1.0 x 10'4n/em?. The de- 
vices under test (DUTs) were SRAM and EEPROM memo- 
ries. Test data transfer was carried out via USB connection 
to a personal computer. The test system allowed for rapid de- 
tection and storage of SEU events. Before radiation, 55 HEX 
(01010101) was written in each storage cell for the DUTs. 
During radiation, the DUTs were read continuously and the 
SEU test could be easily handled by the error detection and 
record system. 


2.0x10° 
1.8x10° 4 —2— 0.18m HM628512C 
s] —O— 0.35m HM62W8512 
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Fig. 4. Cumulative SEUs versus neutron fluence for HM62W8512B 
and HM628512C. 
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Fig. 5. Cumulative SEUs versus neutron fluence for HM628512A. 
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Simulation results predict that SEU for SRAM will oc- 
cur in a reactor neutron environment. The experimental re- 
sults are consistent with the prediction. Figure 4 shows 
single bit upsets induced by neutron irradiation for Hi- 
tachi HM628512C and HM62W8512B. Because of a much 
lower number of SEUs recorded under irradiation for Hitachi 
HM628512A, results for this memory are depicted separately, 
in Fig. 5. The SEU occurs once SRAM is exposed to neutron 
irradiation. In addition, the number of upsets increases lin- 
early with neutron fluence. The SEU cross section is finally 
estimated from the following formula: 


N 
Cae (1) 
where No is total bit number of the SRAM block and is equal 
to 21° bits, N is number of upsets, œ is neutron fluence, and 
N/¢ is equal to the slope of zero-crossing linear fitting. 

The cross section is 3.3 x 10~'4cm*/bit for 0.18 um 
SRAM HM628512C, 4.4 x 10715 cm?/bit for 0.35 um 
SRAM HM62W8512B, and 1.7 x 10716 cm?/bit for 0.5 um 
SRAM HM628512A. As shown in the figure, the total 
SEUs of the HM628512C reached 1.7 x 10° at the end of 
irradiation, reflecting extremely high sensitivity to reactor 
neutrons. The HM628512C shows greater sensitivity than 
HM62W8512B and HM628512A. The reason is that smaller 
feature size sharply increases the upset cross section by low- 
ering the critical charges of memory cells. 

Simulation results predict SEU for EEPROM is impossi- 
ble in reactor neutron environments. SEU tests of EEPROM 
AT28C256 and AT28C040 were carried out in the neutron 
radiation environment. As expected, neutron radiation in- 
duced SEU effects were not significant for EEPROM while 
performing read operations during irradiation, since no error 
bit was observed in lengthy validation tests during irradiation. 


IV. CONCLUSION 


A simulation approach is presented to predict SEU occur- 
rence induced by neutrons in memory devices. The simula- 
tion is adept at processing neutron transport in semiconductor 
materials, obtaining type and energy distribution of secondary 
ions, and calculating LET values. The whole simulation pro- 
cedure provides an approach to get the energy distribution and 
LET spectrum of neutron induced secondary ions, and intro- 
duces an applied method to predict whether SEU will occur 
through comparing the SEU threshold of a memory device 
with the maximum LET in a neutron radiation field, when the 
device is exposed to neutron irradiation with a complex spec- 
trum. The experimental results show three types of SRAM 
devices exhibit high sensitivity to SEU effects, while two 
types of EEPROM devices exhibit no SEU effects. A good 
agreement is obtained between experimental results and sim- 
ulation prediction. For hardness assurance, SEU effects must 
be taken into account when the maximum LET in the radia- 
tion environment exceeds the SEU threshold of the memory 
devices. 
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